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Abstract: 

The recent calculation of the four-loop /^-function in QCD provides further evidence 
that the Banks-Zaks expansion in 16^ — nj is sufficiently well behaved to be useful even for 
n/ ~ 2 light flavours. This expansion inherently predicts "freezing" of the QCD couplant 
at low energies as a perturbative effect. We consider the e~^e~ and Bjorken-sum-rule cases 
as examples. 



1 Introduction 



The old idea that the QCD running couphng "freezes" at low energies has been phe- 
nomenologically successful in a wide variety of contexts. (See Refs. [||, and references 
therein.) Theoretical evidence that "freezing" does occur, and for purely perturbative rea- 
sons, comes from the third-order calculation of Re+e- HI which, when "optimized" with 
respect to renormalization scheme, yields such behaviour Fade approximant methods 
indicate a similar conclusion [Q]. Another approach is the Banks-Zaks (BZ) expansion 
11) §) 0) §) within which "freezing" is natural and ubiquitous. The relevance of the BZ 
expansion to low-energy QCD phenomenology hinges on an extrapolation in the number 
of (massless) flavours, nj, from 16| down to about 2. Our point in this paper is that the 
credibility of this extrapolation has been significantly enhanced by the recent calculation 
of the QCD (3 function to 4 loops 

The BZ expansion ||5|-§ is an expansion in nj^^^ — Uf, where n™* is the critical number 
of flavours at which asymptotic freedom is lost; i.e., where the first /3-function coefficient 
changes sign. In QCD n^''* is 16|. For nj slightly less than 16| the /? function starts out 
negative, but quickly turns positive, and there is thus an infrared fixed point a* close to 
the origin |^. Frovided the couplant a{fi) initially lies between and a* (necessary 
for asymptotic freedom), the couplant will "freeze" to a* as the energy scale ^ tends to 
zero. Thus, in the BZ expansion, "freezing" is natural and universal and will occur for all 
nf < 16|, unless or until the BZ expansion breaks down. 

Ref. § has suggested that the BZ expansion is qualitatively relevant to the real world 
where only 2 or 3 quark flavours are light compared to the QCD A scale. The crucial issue 
is whether the expansion shows reasonable numerical convergence when nj ~ 2. Assuming 
good behaviour, Ref |^] made a prediction for a certain coefficient. As we discuss, this 
prediction is borne out by the recent calculation of the QCD /? function to four-loops |^]. 

The plan of the paper is as follows: Section 2 presents the notation and extracts the 
relevant coefficients. Sections 3 and 4 discuss the BZ expansion for quantities at Q = and 
at general Q, respectively. An important role is played by 7*, the slope of the /? function at 
the fixed point, and we discuss its coefficients (the "universal invariants" of Grunberg [^) to 
fourth order. Sect. 5 contains concluding remarks. Appendix A summarizes the results for 
a general SU{N) colour group. Appendix B discusses the issue of renormalization-scheme 
(RS) invariance and the notion of 'regular' and 'irregular' schemes. The Higgs-decay case, 
which we argue is not a good guide to infrared behaviour, is discussed in Appendix C. 
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2 Notation 



We write the (3 function in the form: 

13 (a) = /i^ = -ha^ (l + ca + C2a^ + csa^ + .••), (1) 

where a = as/ir. The coefficients, in the MS scheme, are [1^, 11, 

/3o = 26 = 11-^^/, 
38 

= She = 102 - yn/, 

2857 5033 325 n 
/32 = 326c2 = —^f + —^)^ (2) 



2 18 ' 54 
/ 149753 „.„_\ / 1078361 6508 



/33 = 1285C3 = —— + 3564C3 - + 



V 6 /V 162 27 

/ 50065 6472 , \ o 1093 o 

+ \ Ca "-f H '^f- 

V 162 81 ^ ) ^ 729 



Here Cs is the Riemann zeta-function (Cs = 1.202056903 . . ., Cs = 1.036927755 . . .). 

For Uf just below 16i, the (3 function has a zero at a* ~ —^j and a* is asymptotically 
proportional to (161 — rif). Its limiting form: 

ao = Jy^^^^"""^) (3) 

serves as the expansion parameter for the BZ expansion Q. Because the constant of 
proportionality is so small, oq remains small (< 0.4) even with rif = 0. To proceed, 
one re- writes all perturbative coefficients, eliminating rij in favour of ag. The first two 
/3-function coefficients, which are RS invariant, become: 

b = — ao, (4) 

c = + -r- (5) 

ao 4 

Within the class of so-called 'regular' schemes [0, |8|, which includes MS, perturbative 
coefficients have a polynomial dependence on nj, and we may write 

Ci = — (ci-i + CifiUo + Ci^ial + ...). (6) 
ao ^ ^ 

The coefficients, in MS, are collected in the table below. 
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C3 



18.5 
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The BZ expansion can be applied to any perturbatively calculable physical quantity 
of the form: 

7^ = a (1 + ria + r2a^ + r^o? + ...). (7) 
In a 'regular' scheme the coefficients rj are polynomials in nj, and hence in oq: 



rifl + fi^iao + rj,2ao + • • • • 



(8) 



Note that a term ?'j,jaQ or oijo^ can be assigned a degree i + j — p, and all terms in any 
formula must have matching degree. 

The prototypical example is the e^e~ ratio: 



hadrons) 



(g) = 3Eg2 (1 + 



(9) 



where, neglecting quark masses, we can write R^+e- {Q) = ^Sg/ (1 + 7^e+e-)i where T^e-t-e- 
has the form (^). [Actually, for 7^e+e- there is a problem in that r2 involves a term 
1.2395(S(7j)^/ (3S(/j^) whose Uf dependence is ambiguous because it depends on the electric 
charges we assign to the additional, ficticious quarks. This arises because Re+e- involves 
not just QCD, but its coupling to electromagnetism. Fortunately, this term seems to make 
only a small numerical contribution. We shall ignore it henceforth.] The coefficients, in 
MS with the renormalization scale ^ equated with Q, are collected in the table below 0]: 



Coefficients in 7^e+e- 



■• [MS(/x = Q)] 
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Another example is the Bjorken sum rule: 

[ dxgf'''"{x,Q'^) = \ 
Jo o 



9V 



(10) 



/o - 3 

(The same QCD corrections, apart from a (Sgj)^/(3Sgj^) term, appear in the Gross 
Llewellyn-Smith sum rule.) The coefficients, from Ref. |jl2| are listed below. 



Coefficients in T^bj [MS(;U = Q)] 



_ n 

12 
107 



-0.917 
13.38 



r2,o 
r2,i 

r2,2 



1385 _ 55 /- 
72 4 



( 107\ ( 2749 I 61 A 

\ — ) + T8^3 - 5C5 

107V /iisA 
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-35.76 
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We mention that the same decomposition of coefficients is needed in the "large-?)" ap- 



proximation |13], which employs the opposite limit ib oo, rather than b = (107/8)ao — > 
as here). 



3 BZ Expansion: Q = 



The fixed-point condition f3{a*) = always has a solution as a power series in oq: 

a* = ao 

A straightforward calculation yields: 



1 + viao + ^200 + ^sflo + • ■ 



(11) 



^^1 = Ci,o + C2-l, 

V2 = (ci,o + 2C2 -l)(ci,o + C2 -l) + C2,0 + C3 _i, 

^3 = cf + 6Ci oC2 -1 + Cl,o(3C2,0 + 4C3 _i -M0C2 _i) 

+ C2 -l(4c2,0 + 5C3 _i) + 5C2 „i + C2,l + €3,0 + C4,. 



(12) 



(Numerically, vi = 1.1366, V2 = 23.27, = 04^-1 — 138.6, in the MS scheme. The poor 
apparent convergence of the a* series need not bother us, since a* is RS dependent.) 

A physical quantity TZ will also have an infrared limit given by a power series in oq. 
One simply takes the perturbative expansion of TZ, Eq. (|^; substitutes a = a* , given by 
(12); and re-expands in powers of qq. This yields: 



TV 



ao 



1 + wiao + W2al + w^aQ + ©(ag) 



(13) 
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where 



,2 



(14) 



v-i + {2v2 + ^1)^1,0 + i'i(2ri,i + 3r2,o) + r-2,i + rs^o- 



These coefficients are RS-scheme independent (see Appendix B) and so their numerical 
values are significant. They should be order-1 numbers, if all is to be well. 

For the e'^e~ case, Ref. Q obtained the value of the first coefficient, wi = 1.22, but 
W2 could only be obtained as —18.25 + c^-i, since 03^1 was then unknown. To quote 
Ref. "For the expansion to be credible one needs C3^_i(MS) to be in the range, say, 
+13 to +21." This prediction is confirmed by the new /3-function result |Q, which yields 
C3 _i = 18.5. Therefore, W2 is quite small, 0.23, giving a respectable series: 



ao 



1 + 1.22ao +0.23ag + 



The next coefficient. 



W3{e+e-) = Ca-I + r3 o(e+e-) - 164.0, 



(15) 



(16) 



would require calculation of both (3 and 7^e+e- to one more order. 
In the Bj-sum-rule case the corresponding result is 



1 + 0.22ao - 1.21a5 + . . . 



(17) 



where the coefficients are also of order unity. It is interesting that the TZ* in this case is 
even smaller than in the e^e~ case. The next coefficient is 



^i'3(Bj) = C4,-l + r3 n(Bj) - 203.7. 



(18) 



4 BZ Expansion: Nonzero Q 

A formulation of the BZ expansion for quantities at a general Q was derived in Ref. |^]. 
We briefiy review the main ingredients. First, we need a suitable form of the boundary 
condition for the /3-function equation. Setting f3{x) = [3 (x) /6, we write 



,in(e 

A 



lim 

5^0 



dx 
s P{x) 



+ C{5) 



(19) 
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The constant of integration C {5) needs to be suitably singular as 5 ^ and we choose 

MS-- 

1"^ dx 



C{5) = P.V.^ 



x'^ (1 + cx) 



= \+c\n6 + c\n\c\ + 0{5). (20) 



Note that Cauchy's principal value (P-V.) is introduced to deal with the pole at x = — 1/c 
when c < 0. This choice amounts to a definition of A, within a given RS. [We use a tilde 



to distinguish it from the older, but still widely used, definition of the A parameter 
The relation is ln(A/A) = (c/&)ln(2 | c| /b). While the two definitions are not dissimilar 
for small Uf, they become infinitely different as n/ ^ 16i. In the BZ-expansion context 
the use of A is much more convenient.] 

As explained in Ref. [^, it is convenient to put the /? function into the form 

^ \ + --TT7^.+H{x). (21) 



^ (x) X 7* (a* — x 

where 7* is 7*/6, with 7* being the slope of the (3 function at the fixed point: 

= -ba* (1 + 2ca* + 3c2a*^ + 4c3a*^ + ...). (22) 
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dx 

As discussed below, 7* can be obtained as a series in oq- The remainder function H{x) 
can be expanded as a power series, Hq + Hix + . . ., whose coefficients are of order ao- 

One now inserts ( |2l| ) into (|l9|) and performs the integration. One can then eliminate 
a and a* in favour of IZ and TZ* . In fact, since the result must be RS invariant, one can 
— without loss of generality — short-cut this step by utilizing the "effective-charge" RS 
in which a = TZ. This leads to the formula |^]: 

= ;^ + ^ln 1--^ +cln(|c|7e) + V \ , ■ (23) 



The last term, involving the H^^'^^ coefficients (of the effective-charge scheme), is only 
relevant in fourth order and beyond. Thus, for the first three orders the equation takes 
the same form, just with the parameters 7* and TZ* approximated to the appropriate order. 



On the left-hand side, pi is the RS invariant [14| 

(I) „(£-), (24) 

where Agfr is a characteristic scale specific to the particular physical quantity TZ. It is 
related to the A parameter of some reference scheme (eg. MS) by an exactly calculable 



6 



factor exp(ri/6) involving the ri coefficient in that scheme, evaluated ai fi = Q. (We 
caution that pi cannot be split into 0(1) and ©(ao) pieces in a RS-invariant way [^.) 

Numerically inverting Eq. (|23| ) provides 7^ as a function of Q. The resulting 1Z{Q) 
naturally agrees with ordinary perturbation theory to the corresponding order at large 
Q, but freezes to the value 7^*(") as Q ^ 0. The BZ series expansion for TZ* was dis- 
cussed in the previous section. The BZ expansion for 7* is obtained straightforwardly by 
substituting the expansion of a* (Eqs. (Ill]) and (p^) into ([2^). This gives: 



7 = ao 



1 + giao + 52ao^ + 53^0^ + • • • 



(25) 



where, 



91 = ci,o, 

92 = ci/ - C2-1' - Oi-i, (26) 
53 = ci,o^ - 4C2 -1^ - 5C1,0C2 -1^ - 4ci,0C3 -1 

- 2C2,-.lC2,0 - 6C2 -IC3 _i - C3,o - 2C4 _i. 



It is noteworthy that certain terms of degree n are absent in gn'- 91 does not contain C2,-i; 
52 does not contain C2,o or C2,-iCi^o; and 53 does not contain C2,i or C2,oCi,o or C2,-iCi^o^- 

The significance of 7* = 67* is that it is the 'critical exponent' governing how TZ 
approaches 7^* as Q — > 0; asymptotically, TZ — IZ* (x. Q"' . As pointed out by Grunberg, 
the gn coefficients are RS invariants, and are universal, in the sense that they are not 
specific to some particular physical quantity TZ. (See Appendix B for discussion of some 
subtleties.) 

The new /?- function result Q enables us to determine the numerical value of the second 
invariant; g2 = —8.99. (The exact expression, for general A^, is given in Appendix A.) 
Hence the 7* series is: 

J* = ao [1 + 4.75ao - 8.99ao^ + . . .] (27) 

Clearly, the 7* series is not as well behaved as the TZ* series that we saw earlier. In Fig. 
1(a) we show 7* (= 67*) as a function of nj. (Fig. 1(b) shows shows the same quantity 
normalized by l/(6ao).) The lower and upper solid curves are the first- and second-order 
results, respectively, while the middle solid curve is the third-order result. The dashed 
curve represents the third-order result re-cast as a Fade approximant: 
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For comparison we also give, as the dotted curve, the prediction arising from an optimized- 
perturbation-theory analysis of the e~^e~ case jl[ (see comments in Appendix B). The 
reasonable agreement between the last three curves gives us some confidence that the 
extrapolation to low is qualitatively valid, even if the quantitative precision is not 
good. 

The next coefficient is 

53 = 269.44 - 2c4-i. (29) 
We therefore predict that C4,_i(MS) will turn out to be somewhere around 135 it 10. 



Assuming that the W3 coefficients in ( |lq) and (p^ are modest, we can expect that 
^3,o(e+e-) 29 ± 10 and r3,o{Bj) ?a 68 ± 10. 

Knowing 7* and TZ* , we can numerically solve Eq. (]23|) to obtain 'JZ{Q) as a function 
of Q/Acff . The result, for the e~^e~ case in third order, is shown in Fig. 2 for Uf = 14, 10, 6 
and 2. (It may be directly compared with similar figures for first and second order in Ref. 
Q.) Corresponding results in the Bjorken-sum-rule case are shown in Fig. 3. Note that 
the Aefj's defining the units of Q in Figs. 2 and 3 are not the same: however, they are 
easily converted to a common A using Eq. (p^. 

5 Concluding Remarks 

The result of van Ritbergen et a/.'s |^] calculation of the /? function to four- loops sheds 
much light on the BZ expansion. It supports the idea that the expansion is relevant to the 
phenomenologically interesting case of only two light quark fiavors. That, in turn, implies 
perturbatively explicable "freezing" of the QCD running coupling constant. 



We mention again the very interesting work on the "large-6" approximation (see |13] 
and references therein) which draws on large-nj results. This is the opposite approximation 
to ours. It extrapolates upwards from nj = —00 {minus infinity if the theory is to be 
asymptotically free) towards nj ~ 2, whereas we are extrapolating down from nj = 16i. 
We think both approximations are useful; neither is very precise, but both seem to be 
qualitatively valid, and offer a great deal of insight into QCD. Our preliminary studies 
indicate that the large-6 approximation also predicts "freezing," and we hope to report on 
this shortly. 

Acknowledgements: We thank Jifi Chyla for correspondence, particularly regarding 
the issues in Appendix B. This work was supported in part by the U.S. Department of 
Energy under Grant No. DE-FG05-92ER40717. 
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Appendix A: SU(A/^) generalization 

The critical number of flavours is: 



11 



nf^ = ^N. (30) 

The /3-function coefficients in MS are 

/3o = 26 = ^(lliV-2n/) 
/3i = 86c = ^(34iV3-n/(13iV2_3)) 

/32 = 326cp = — 3-^f5714iV^ + n.(-3418iV^ + 561iV2 + 27)+n2.(224A^^-66iV 
108A'^ V 

/33 = 1286c^^ = Y^^^^(601892iV^- 25920iV^ + C3(95047V'^ + 6842887V^) 



485513iV*^ + 58583iV^ - 21069iV^ - 5589 
+ C3(-4320iV^ - 1183687V^ + 9504iV2) 

+ n) [69232iV^ - 198167V^ - 224287V + C3(181447V^ - 13824iV3 + 52704iV) 
+ n)[10407V^ - 616iV^]) (31) 



The BZ expansion parameter becomes 

3(25iV2 _ II 

The invariant coefficients are 



«o= ,.o...2 ..M ''-^f) (32) 



(25jV2 - 11) 

= 77 «o, (33) 



16 

ao ' 8iV 



1 (13iV2_3) 

+ ... • (34) 



In MS 



(-1402A^^ + 242N'^ + 33) 

"^''-^ ~ 48iV(25iV2 - 11) ' ^ ' 

(318iV^ + 55jV2-9) 
C2,o = ^^^^^^ , (36) 

(25iV2 - ll)(112iV2 - 33) 

'^•^ - 3072iV ^^^^ 

(14731iV6 - 300477V^ - 58839iV2 _ 2277) 
^^'"^ ~ 1152A^2(25iV2 _ 11) 

11 (25jV4 - i8iv2 + 77) 
+ Y (25iV2 - 11) (^^^ 



The Grunberg invariants, the coefficients in the expansion of 7* are: 
^ (13iV'-3) 

(366782iV^ - 865400A^*^ + 1599316A^^ - 571516iV2 _ 3993) 
~ 768Af2(25iV2 -11)2 

_ll(25iV^-18iV2 + 77) 

8 (25iV2_ii) ^^"^ 

Appendix B: Regular and Irregular Schemes 

In this appendix we discuss renormaHzation-scheme invariance of the BZ expansion results. 
Since the expansion parameter 00 = 3I1 (16^ — rif) is an RS-invariant pure number, one 
expects the coefficients in the BZ expansion of a physical quantity TZ to be RS-invariant, 
but it is important to have confirmation. We first consider the fixed point results at Q = 
(which are independent of the A parameter). 

In the text, as in Ref |^], we limited the discussion to so-called 'regular schemes' in 
which the coefficients in TZ and /3(a) have a polynomial dependence on nj. Such schemes 
are natural in diagrammatic terms, since each fermion loop gives an Uf factor. They are 
convenient for our purposes, since the BZ-expansion coefficients are easily extracted from 
calculations made in those schemes. However, we emphasize that "irregular" schemes are 
not necessarily "bad." They will also lead to the same BZ-expansion results, but the 
extraction of the BZ coefficients from calculations in those schemes will be less straight- 
forward. An analogous situation arises, for example, with gauge invariance; a certain class 
of gauges may be convenient for some purposes, but inconvenient for others. 

It is easy to show the invariance of the BZ coefficients within the class of 'regular' 
schemes. This was done in Ref. |8| by considering the RS-invariants P2,P3, ■ ■ ■ |14] that 



are invariant combinations of TZ and /3(a) coefficients, and expanding them in powers of 
uq. a 'low-brow' version of the proof is also instructive. Consider a change of RS, a a', 
where 

a = a (^1 + uia + U2C? -|- . . . j . (41) 
The coefficients in the expansion of TZ change to 

A = rx-ux, 

r'2 = r2-U2- 2uiri + 1u\ , (42) 
etc.. The /3 function transforms to f3'{a') = (da' /da) (3 {a), whose coefficients are: 
c' = c 



10 



C2 = C2+U2 -uf - UlC, 

C3 = + 2u3 + cul - 2c2Ui - 6uiU2 + 4:ul, 



(43) 



etc.. If both the primed and unprimed scheme are regular, then the Cj, c[ and rj, r'- 
coefficients are expandable as in Eqs. (|6|), (^), and the Ui coefficients in the scheme 
transformation can be expanded as: 

^^1 = ^1,0 + ^^1,100, 

U2 = U2fi + ^2,1^0 + ^2,200 , (44) 

etc.. It is then straightforward to prove the invariance of the combinations appearing in 
the BZ expansion of TZ and 7*. For instance, from the relations 

'^'1,0 = n,o - ""1,0 , 
C2-I = C2-I + ^1,0 , 

c'3-1 = C3 _i - -Ui o - 2c2,-iMi,o, (45) 

one sees that r'lQ + C2_i = ri^ + C2,_i, and c'^-i + (c2^_i)^ = c^-i + (c2,-i)^, showing 
that these combinations are RS invariant. Extending this procedure one can prove the 
invariance of the higher-order Wi and gi coefficients. 

Any scheme related to MS(/i = Q) by a transformation (|4l| ) with Uj's expandable 



as in (44) is a 'regular' scheme. [It is noteworthy that in a general, regular scheme the 
coefficients C2.2,C3^3, . . . are non-zero, while in a more restrictive class of 'strictly regular' 
schemes, of which MS(/i = Q) is an example, these coefficients vanish. This distinction is 



unimportant for the BZ expansion, but matters for the large-6 approximation [13|.] The 
'effective charge' (or 'FAC') scheme, in which a = TZ for some specific physical quantity, is 
a 'regular' (but not 'strictly regular') scheme. However, it is easy to construct RS's that 
are 'irregular' simply by considering a transformation in which the Uj depend on nj in 
a non-polynomial fashion, so that the coefficients are no longer expandable in positive 
powers of gq. There is nothing intrinsically 'bad' about such schemes. They can arise quite 
naturally. For example, the 't Hooft scheme, in which C2 = C3 = . . . = 0, is 'irregular' 
p!3[ |. [Recall that, in 'regular' schemes, C3-1 + (c2,-i)^ is invariant and does not vanish.] 
The principle-of-minimal-sensitivity (PMS) scheme for any given physical quantity is also 
'irregular.' In both these cases the rj coefficients have 1/ao pieces. 

In a 'regular' scheme, obtaining the BZ expansion to nth order requires terms of order 
n + 1 in the /? function, and of order n in the physical quantity TZ. In an 'irregular' scheme 
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the same information is distributed among higher-order coefficients as weh. Starting from 
an 'irregular' scheme, one would require some knowledge of perturbative coefficients to 
higher orders; maybe to all orders. However, provided one carefully kept all terms that 
could contribute to a given order in oq, one would obtain the same BZ expansion. 

Finally, we discuss the finite-Q case. Following Ref. P], we have formulated the result 



as Eq. (23), which is to be solved numerically to obtain 7^ as a function of Q. This 
formula involves the RS-invariant quantities pi„ 7*, TZ* , c, and h\^^\ The A parameter 
and Q appear only in pi. The h\'^^^ coefficients (relevant only in 4th order and beyond) 
are directly related to the p2,P3,--- invariants of Ref. ||l^ which can be conveniently 
redefined (hence the tilde) to coincide with the /9-function coefficients of the 'effective 
charge' scheme |16|. 

The RS invariance of the BZ expansion of 7* is verifiable by the procedure dis- 
cussed above. It is expected because of the well-known result [ |17| | that the slope of the 
P function at a fixed point is an invariant. However, there is an important caveat to the 
last statement which necessitates some further discussion. The quoted result follows 
by differentiating the /3-function transformation, (3'{a') = {da' /da)f3{a), to give 

The first term vanishes at the fixed point — provided that neither da /da' nor d'^a' /da^ is 



singular there |17|. Chyla |18| ] has pointed out that, in general, it can be quite natural for 
those factors to be singular (and arbitrary scheme transformation can of course make fixed 
points appear and disappear!). In general, then, the critical exponent 7* in TZ — TZ* oc Q"^ 
as Q — > is not the same as dfi / da\a=a* ■ However, these two quantities will coincide in 
a large class of schemes. It also seems safe to assume that they coincide in the context of 
the BZ expansion, where there is necessarily a fixed point at a* = a^ + . . .. 

[We find that the 3rd-order PMS scheme, though 'irregular,' also yields 7* = d/3/da\a=a* 
This requires a detailed analysis of the optimization equations |Q, |^ as Q ^ 0. The PMS 
results (in the e~^e~ case) for 7* are shown as the dotted curve in Fig. 1.] 



Appendix C: Higgs decay 

In this appendix we discuss the case of Higgs-boson decay into hadrons. This seems to be 
a much more problematic than the cases discussed earlier. It could be viewed as confiicting 
with the general picture we have presented. We shall argue, though, that the crucial role 
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of quark masses in this quantity makes it unsuitable as a guide to the infrared behaviour 
of perturbative QCD. First, let us discuss the numbers. 

The hadronic decay width of the Higgs has the form Th = 4^ "^5^(0) with 

r(a) = 1 + Tia + . . ., and where rriq is the running quark mass evaluated at some scale. 
For Th there is a factorization-scheme ambiguity (how much of the radiative corrections 
should be absorbed into m^, and much should be left in the explicit series r(a)?). However, 
one can define the quantity |19] 



T^Hi: 



ggs 



2 dlnMjj 



(47) 



which is free of this factorization-scheme ambiguity, and is a physical quantity of the same 



form TZ = a(l + ria + . . .) considered earlier. The coefficients, from Ref. [19| are collected 
in the table below. 



Coefficients in T^Higgs- [MS(;U = Q)] 



23 

12 



( WL\ ( 11 

V 8 M 6 



1.917 
24.52 



r2,o 
^2,1 

^^2,2 



503 _ 55 A 

18 4 ^3 
107 A ( 2935 _ 7 A _ 7i£ 

8 y \^ 144 4''3 3 



107\ ( 275 /• 



12 



-44.47 
200.5 
321.1 



^^3,0 
^3,1 
^^3,2 



22631621 
82944 

mi: 



432 
208411 



3456 

107 A ^ ( 694303 _ 
8 y V 6912 
107 ( 985 5 
8 y VIO8 2 



C3 + -^Q5 + 

20755 a _ 355 a _ 137„2 
S3 4« "55 1Q2 " 



288 
119 



48 

4 C3 + 
C3 - 24^ 



317^2 
96 " 



-166.6 
-2162 

6901 

3808 



In Ref. [20 1 it was observed that, at third order in the effective-charge (or 'FAC') 
scheme, there is a fixed point with T^Jjiggs ~ a* ~ 0.15. The authors viewed this as 
probably spurious. Indeed, it is odd that it is only about half the size of the frozen 
couplant found in the e^e~ case, and so is far from the leading-order BZ expectation that 
a* ~ cq. At 4th order Ref. |1£] finds that this fixed point is no longer present. We have 



checked that the situation is much the same in optimized perturbation theory [14, 21 1 
The (5 = BZ series in this case is: 



^Higgs = oo(l + 3.05ao + 7.67ao + . . .) 



(48) 



whose coefficients are considerably larger than in the e~^e or Bjorken-sum-rule cases (Eqs. 
(15), (p^)). For a low number of flavours the "corrections" are as big as the leading term, 
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and both the same sign. (The next coefficient in the expansion is C4^_i — 109 and is still 
unknown; our estimate (see Eq. ( [29| ) suggests that it is around 26 ± 10.) One could 
conclude that this is perhaps a case where the BZ expansion breaks down before nj = 2, 
so that maybe there is no freezing in this case. 

We believe, however, the problem is that this quantity is just not a useful indicator of 
massless QCD's infrared behaviour, because of the way that quark masses are involved. 
For exactly massless quarks the hadronic Higgs-decay rate is zero, because the Higgs- 
quark coupling is proportional to quark mass. (The calculations neglect quark masses in 
the radiative corrections but not, of course, in the overall coupling factor.) If we keep the 
quark masses finite when we consider, theoretically, the limit Mh — s- 0, we will trivially get 
zero as soon as the decay becomes kinematically forbidden. To avoid this we would need 
to consider a limit in which ruq tends to zero at least as fast as Mh', say niq oc (Mh)'^ with 
K > 1. But then T^Higgs is not of the form a(l + ria + ...), and depends on k, making it 
ill-defined. These issues do not arise for 7^e+e- ^Bj) which are meaningful for massless 
quarks. 
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Fig. 1. (a) The critical exponent 7* in first, second, and third orders of the BZ expansion 
(lower, upper, and middle solid curves). A Pade approximant form of the third-order result 
is shown as the dashed curve. The dots represent the result of an optimized-perturbation- 
theory analysis [1]. (b) The same, normalized by l/(6ao); i.e., 7*/ao. 
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Fig. 2. 7?.e+e- as a function of Q/Aes to third order in the BZ expansion 
forn/ = 14, 10,6,2. 
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Fig. 3. TZb} as a function of Q/Ags to third order in the BZ expansion 
forn/ = 14, 10,6,2. 
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